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The structural, electronic, and magnetic properties of the superconducting NdFeAsOi-^F^ phases 
(Tc — 43 K for x — 0.15) have been investigated experimentally by high-resolution synchrotron x-ray 
powder diffraction, magnetization and 75 As NMR measurements. Density-functional calculations 
were performed to calculate and analyze the electric field gradient and the density of states. Com- 
pared to LaFcAsO family, the NdFeAsO family shows a contraction of the lattice parameters with 
shorter rare-earth (RE)-As bond distances, an increased thickness of the As-Fe2-As layer, and less 
distorted Fe-As4 tetrahedra. The 75 As quadrupole frequencies are enhanced with respect to the La 
analogs. This is due to a more prolate As 4p electron distribution mainly caused by the reduced 
lattice parameters and not by the presence of Nd 4/ electrons. A non-negligible hyperfine cou- 
pling between the 75 As nuclei and the Nd 4/ states indicates a weak coupling between the REO 
and FeAs layer and possibly opens the channel for a Ruderman-Kittel-Kasuya-Yosida (RKKY)-type 
interaction between localized Nd 4/ moments mediated by itinerant Fe 3d and/or Nd 5d states. 

PACS numbers: 74.70.-b, 76.60.-k 



I. INTRODUCTION 

The family of rare-earth (RE) quaternary oxypnictides 
with the general formula REFeAsO (RE = La, Cc, Pr, 
Nd, Sm, Gd, Tb, Dy) (Ref. [0) is currently generat- 
ing much excitement in the condensed matter commu- 
nity. When electron doped (e.g., partial replacement 
of 2 ~ with F~ or oxygen deficiency), these materi- 
als are reported to display superconductivity with a su- 
perconducting transition temperature, Tc, as high as 
55 K,2 & i £&L £ & 12 In electron-doped LaFeAsO the su- 
perconductivity (Tc ~ 26 K) appears in the close vicin- 
ity to a spin-density wave ground state ffi 11 ' 12 implying 
that antifcrromagnctic (AF) correlations may be impor- 
tant in the promotion of the superconducting ground 
state. La 3+ has no unpaired electrons, and the mag- 
netic response can be attributed solely to the Fe 3d 
electronic states that dominate the region around the 
Fermi level. The relation between AF fluctuations and 
superconductivity has been proposed theoreticall y 13 i 14 
and addressed experimentally by 19 F, 75 As, and 139 La 
NM R 15 ' 16 ! 17 . Upon F doping, a pseudo-gap behavior has 
been found in the temperature dependence of the 19 F 
(Ref. GJ) and 75 As (Ref.© 

spin-lattice relaxation rates. 
Substitution of La 3+ with other rare-earth elements 
such as Nd 3+ (free-ion magnetic moment /i ~ 3.6 /xg) in- 
troduces additional spin degrees of freedom in the REO 
layer due to the unpaired RE 4/ electrons. Interest- 
ingly, the NdFeAsOi-^Fa; family displays superconduc- 
tivity with enhanced Tc's of ~ 50 K compared to the La- 
based materialsi2i£ This raises important questions of if 
and how the rare-earth magnetic moment interacts with 
the conducting FeAs layer and whether these 4/ moments 
or lattice effects (chemical pressure) have a stronger influ- 



ence on the electronic structure and the related magnetic 
and superconducting properties in the REFeAsO family. 

To address these questions, we have studied the struc- 
tural, magnetic, and superconducting properties of the 
NdFeAsOi-zF^ family by means of high-resolution syn- 
chrotron x-ray powder diffraction, magnetization, and 
75 As NMR measurements. Our experimental results are 
supported by electronic structure calculations. 



II. METHODS 

Polycrystallinc NdFeAsO and NdFeAsO0.85F0.15 sam- 
ples were prepared by a standard two-step solid-state re- 
action method. The magnetic susceptibility was mea- 
sured using a Quantum Design magnetic property mea- 
surement system. Synchrotron powder diffraction data 
sets, collected on the ID31 high-resolution diffractome- 
tcr at the European Synchrotron Radiation Facility in 
Grenoble, France, were binned with a 0.002° stcpsize 
(A = 0.309952 A). The GSAS suite of programs and 
EXPGUI graphical user interface were used for Rietveld 
fitting, 1 ^ 75 As (I = 3/2) NMR frequency-swept spec- 
tra were measured in 8.9 T magnetic field with a selective 
solid echo pulse sequence (pulse length of 7r/2 = 22 /zs 
and delay of r = 50 (is). The reference frequency of 
i/( 75 As) = 65.096 MHz was determined from a NaAsFg 
standard. 

The band-structure calculations were performed us- 
ing the full-potential local-orbital minimum basis code 
FPLO (version 5.00-19) (Ref. S3) within the local den- 
sity approximation(LDA). In the scalar relativistic calcu- 
lations the exchange and correlation potential of Pcrdew 
and Wang2i was employed. To treat the Nd 4/ state 
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adequately LSDA+(7 in the atomic limit was employed 
with [/ Nd = [7,9] eV and J = 0.7 eV. As basis set 
Nd (4p4d4f5s5p/6s6p5d+7s7p), Fe (3s3p/4s4p3d+5s5p), 
As (3p/4s4p3d+5s5p) and O (2s2p3d+3s3p) were chosen 
for semicore/valence+polarization states. The high lying 
states improve the basis which is especially important for 
the calculation of the electric field gradient (EFG) ten- 
sor with components Vij = dV/dxidxj. The lower lying 
states were treated fully relativistic as core states. A well 
converged k mesh of 252 k points was used in the irre- 
ducible part of the Brillouin zone. In order to investigate 
the influence of F substitution on the O site, the virtual 
crystal approximation (VCA) was applied. 



III. RESULTS 
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High-resolution synchrotron x-ray powder diffraction 
profiles of NdFcAsO at room temperature confirmed that 
its structure is tetragonal — space group PA/nmm — in 
analogy to other REFeAsO systems. Rietveld analy- 
sis of the synchrotron X-ray powder diffraction data 
(Fig. QJ resulted in lattice parameters a = 3.96629(1) A 
and c = 8.59886(6) A with Fe-As and Fe-Fe distances 
of 2.4007(4) A and 2.80459(1) A, respectively and an 
As-Fe-As bond angle of 111.39(3)° [\ 2 = 2.4, wRp = 
13.2%, Rp = 9.2%, Nd: z = 0.13887(6), and As: z = 
0.65735(9)] 32 Compared to LaFeAsOfii the Nd 3+ con- 
taining oxypnictide shows a contraction of the lattice con- 
stants with shorter RE-As bond distances, an increased 
thickness of the As-Fe2-As layer and less distorted Fe- 
AS4 tetrahedra. Temperature-dependent synchrotron x- 
ray powder diffraction measurements collected between 
room temperature and 100 K revealed the presence of a 
structural phase transition to an orthorhombic (O) su- 
perstructure (ao ~ bo — V2ax and cq = or\ space 
group Cmma) in analogy with the results reported in 
Ref. [23]. All the Bragg reflections containing the tetrag- 
onal Miller indexes h ^ and k ^ start broadening at 
160 K and finally split only at 140 K (inset to Fig. [J). 

The F-doped compound, with nominal composition 
NdFcAsO0.85F0.15, also adopts tetragonal symmetry but 
with reduced room temperature lattice parameters [a = 
3.96251(1) A and c = 8.57785(6) A)] and Fe-As and Fe- 
Fe distances of 2.4025(6) A and 2.80192(1) A, respec- 
tively (x 2 = 2.4, wRp = 12.84%, Rp = 10.31%, Nd: 
z = 0.14211(6), and As: z = 0.6544(1)). We notice that 
despite the smaller ionic radius of F~ compared to O 2- , 
the RE-0 distance increases, resulting in an increased 
thickness of the Nd-CVNd layer and more regular Nd4- 
O tetrahedra. 

The zero-field-cooled (ZFC) molar magnetic suscep- 
tibility, x, of NdFeAsO measured in a 10 kOe applied 
field shows a substantial Curie- Weiss-type contribution 
[Fig. [2^] that completely dominates \ at low tempera- 
tures. Analysis of the susceptibility data gives the Curie 
constant C = 1.43(2) emuK/mol and the Curie-Weiss 
temperature 9 = -14(1) K [inset to Fig.Hk)]. The Curie 



FIG. 1: (Color online) Measured (O) an d calculated (-) syn- 
chrotron x-ray powder diffraction (A = 0.309952 A, room tem- 
perature) profiles for NdFeAsO. The lower solid line shows the 
difference profile, and the tick marks show the reflection po- 
sitions (top: NdFeAsO, bottom: Nd203 impurity). The inset 
illustrates the T-O structural phase transition. 

constant corresponds to an effective magnetic moment 
fJ-eff = 3.31(4) (j,b in agreement with the expected value 
for localized Nd 4/ moments, while the negative Curie- 
Weiss temperature implies AF interactions between the 
Nd 4/ moments. We also note here that the concentra- 
tion of Nd2 03 impurities is so small that the values of 
C and 9 do not change after molar susceptibility is cor- 
rected for the impurity contribution. The temperature 
dependence of the zero-field-cooled magnetization mea- 
sured on NdFcAsO0.85F0.15 [Fig. HJd] shows a clear onset 
of bulk superconductivity at Tc = 43(1) K with shielding 
fraction of 40%. 

Representative powder 75 As NMR spectra of NdFeAsO 
and NdFeAsO0.85F0.15 arc shown in Fig. [3J The room- 
temperature spectra for both samples display a typical 
quadrupole powder line shape. For undoped NdFeAsO 
the shape of the central —1/2 <-> 1/2 transition and the 
sharp singularities for the satellite transitions strongly 
indicate that the asymmetry parameter, defined by the 
components of the EFG tensor as 77 = (V xx — V yy ) /V zz , 
is T] — 0. V zz defines the quadrupole frequency 

_ 3 eV zz Q 

VQ - 21(21 -iy (1) 

with Q as the quadrupole moment of 75 As. We recall that 
the 75 As nuclei have four nearest-neighbor Fe atoms (in- 
set to Fig. [3|) and in the tetragonal phase reside on the 2c 
(1/4, 1/4, z) position, which is fourfold axially symmetric. 
77 = is therefore fully consistent with the As site symme- 
try. The 75 As powder NMR line shape simulation (inset 
to Fig.[5]) then gives vq = 11.8(1) MHz. We note that v Q 
is significantly enhanced compared to vq = 8.7(4) MHz 
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FIG. 2: (a) Temperature dependence of the ZFC susceptibil- 
ity measured in 10 kOe applied field for NdFeAsO powder. 
Inset: temperature dependence of the inverse ZFC suscepti- 
bility for NdFeAsO powder, (b) Temperature dependence of 
the ZFC magnetization for NdFeAsO0.85F0.15 with the super- 
conducting transition temperature of 43 K and a shielding 
fraction of 40.6%. The magnetic field was set to 20 Oe. 



(Ref. l24h for LaFeAsO and it is even larger than vq = 
11.00(5) MHz reported for LaFeAsOo.gFo.i^ This is pre- 
sumably a result of slightly different charge distributions 
around the As site reflecting the small structural varia- 
tions in the As local environment as it will be discussed 
below. The 75 As NMR spectra of NdFeAsO0.s5F0.15 are 
considerably broader compared to the parent compound, 
reflecting local disorder introduced by the F doping. This 
prevents a precise determination of 77. However, the As 
NMR spectra can be simulated using a Gaussian distri- 
bution of vq with a mean value of vq = 12.9(2) MHz and 
a full width at half height of the Gaussian distribution 
Av Q = 1.8(2) MHz. 

For the undoped compound, the measured EFG is 
obtained by applying the 75 As quadrupole moment 2 ^ 
Q = (0.314 ± 0.006) b in Eq. © and yields \V™ P \ = 
(3.11 ± 0.09) • 10 21 V/m 2 . Using the experimental room 
temperature lattice parameters and atomic positions and 
£/ Nd = 8 cV we obtain a good agreement for the calcu- 
lated EFG: V°f c = -3.39 ■ 10 21 V/m 2 . Varying C/ Nd by 
±1 eV (within the physically reasonable range) changes 
the EFG only by qpO.01 . 10 21 V/m 2 , which is well below 
the experimental error bars. Like for the Fe magnetic 
moment, which in the calculations shows a strong depen- 
dence on the As z position , 26 ! 27 we observe also for the 
EFG a strong As z dependence (inset to Fig. 2]). The cal- 
culated EFG (V°f c = -2.92 ■ 10 21 V/m 2 ) is even closer 
to the measured one when As is shifted along the neg- 
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FIG. 3: Representative 75 As NMR spectra measured in (a) 
the parent NdFeAsO and (b) the F-doped NdFeAsO0.85F0.15 
powder samples. Inset: structure around the As ion show- 
ing its room-temperature distances to nearest-neighboring Fe 
(black spheres) and Nd ions (white spheres). 

ative z direction to z = 0.6515. Here, the energy has a 
minimum (in inset to Fig. [4] marked by an arrow) and 
the structure has a shorter Fe-As distance of 2.3729 A. 

Analogous to the experiments, where the quadrupole 
frequency for 75 As is larger for the Nd than for the La 
compound, the magnitude of the calculated \V ZZ | is larger 
for the Nd compound V zz = —3.39 • 10 21 V/m 2 than for 
the La compound 2 ^ V zz = -3.21 • 10 21 V/m 2 but less 
drastically. Analysis of the calculated EFG (Ref. [29t ) 
shows that the largest contribution to the EFG mainly 
originates from As Ap electrons. In both compounds the 
anisotropy count A p = ^(n x + n y ) — n z for the As Ap 
electrons is negative, which corresponds to a prolate p 
electron distribution and is in agreement with the nega- 
tive V zz . For the La compound A p is less negative than 
for the Nd compound: A p = -0.018 and A p = -0.062, 
respectively, which means that the Ap orbitals are less 
isotropic for the Nd compound and hence the EFG is 
enhanced. The observed difference between the La and 
Nd compounds could arise from cither the additional 4/ 
electrons due to screening effects or from the smaller lat- 
tice parameters in the Nd compound. To separate the 
influence of the change in the lattice geometry from the 
change in the electronic configuration, we compare the 
real EFGs with fictitious EFGs, which are obtained by 
exchanging the lattice parameters for both compounds 
(see Table I). 

The \ V ZZ \ value increases with lattice compression both 
for the LaFeAsO and NdFeAsO electronic configurations 
(Table I from the left to the right). Next, we checked 
the role of the RE Af electrons. In this case we find the 
opposite trend; i.e., for given lattice parameters \V ZZ \ is 
smaller for the Nd than for the La compound (Table I 
from top to bottom). One can therefore conclude that 
I V zz I increases with lattice contraction and decreases with 
increasing number of RE Af electrons. 
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TABLE I: Calculated V zz for As for both compounds for both 
lattice parameters. In the calculations we used the follow- 
ing room-temperature lattice parameters: set 1 [LaFeAsO 
(Ref. Gil)]: a=4.03007(9) A, c=8.7368(2) A , As: z = 
0.6507(4), and La: z = 0.1418(3); set 2 (NdFeAsO, this work): 
a = 3.96629(1) A and c = 8.59886(6) A ; As: z = 0.65735(9) 
and Nd: z = 0.13887(6). LaFeAsO with set 1 and NdFeAsO 
with set 2 are the real compounds; the other combinations 
are fictitious. 
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Now, we focus on the effect of electron doping on 
the EFG in NdFeAsO. The EFGs of the doped com- 
pounds were calculated with the virtual crystal approx- 
imation. The validity of the VCA was confirmed by su- 
percell calculations.— First, we consider solely the ef- 
fect of electron doping. Therefore, we keep the struc- 
tural parameters fixed for different levels of doping. In 
Fig-Stwo such VCA curves are shown. When the exper- 
imentally determined As z = 0.65735(9) position (Nd- 
FeAsO at room temperature) is used, the calculated and 
measured EFG values for 15% doping agree very well. 
The VCA curve with the optimized As z = 0.6515 posi- 
tion deviates from the experimental curve as it system- 
atically predicts a smaller \V ZZ \. Next, we investigate 
the structural change on top of the doping by calcu- 
lating the EFG within VCA for the room-temperature 
structural parameters of NdFeAsO0.85F0.15- We notice 
that the use of these parameters slightly reduces \V ZZ \ to 
V zz = -3.17 ■ 10 21 V/m 2 . Please note the black circle 
marked by an arrow in Fig. 2] We stress at this point, 
that the effect of electron doping on the EFG is much 
smaller than the influence of the As z position as can be 
clearly seen by comparing Fig. [4] with the inset of Fig. [4] 
Our calculations predict a small decrease in \V ZZ \ upon 
electron doping for NdFeAsO0.85F0.15, although experi- 
mentally a slight increase is observed. This is also the 
case for the La compound where the difference between 
the slopes of the calculated and the experimental EFGs 
upon doping 2 ^ is even more pronounced. Further studies 
are required to investigate if this is due to the quality of 
the samples or due to intrinsic changes in the electronic 
structure. 

In undoped NdFeAsO the 75 As NMR line shape starts 
to change below T5 ~ 160 K (Fig. [3]). Its evolution can 
be accounted for by assuming that r\ suddenly starts to 
increase below T5 and reaches 77 = 0.10(2) at 140 K. 
The increase in 77 is a strong indication that the 75 As 
site symmetry has lowered due to the structural phase 
transition to a low-temperature orthorhombic structure 
(space group Cmma)^ similarly as found in other RE- 
OFeAs systems. We note that in the orthorhombic phase 
75 As is on the Ag (0, 1/4, z) position, allowing r) to differ 
from zero. The band-structure calculations show that the 
influence of the orthorhombic distortion on the EFG is 




1.2 -0.1 0.1 0.2 

As z displacement [A] 



l_| , 1 , 1 , 1 , 1 

0.05 0.1 0.15 0.2 

doped electrons 

FIG. 4: (Color online) Calculated V zz obtained from the vir- 
tual crystal approximation using the experimental As z — 
0.65735(9) position (empty squares), the optimized As z — 
0.6515 (filled diamonds) position, and the NdFeAsO0.85F0.15 
structure at room temperature (black circle marked by an ar- 
row) . The measured EFGs for the pure and the 15% F-doped 
compound are shown by error bars. Inset: dependence of 
75 As V zz on the As z position. The optimized As z position 
is marked by an arrow. The experimental V zz is represented 
by the shaded bar. 



negligible - much less than the experimental error bars. 
This is understandable given the small size of the or- 
thorhombic splitting. The anisotropy parameter is more 
sensitive to the distortion and increases to 77 = 0.04 (com- 
pare to the experimental value of 77 = 0.10(2)). Below 
T)v = 140 K the high-temperature 75 As resonance in Nd- 
FeAsO starts to rapidly lose in intensity on the account 
of an extremely broad low-temperature signal (inset to 
Fig. [5]). The dramatic broadening of the 75 As resonance is 
a manifestation of AF magnetic order of Fe moments oc- 
curring at T/y = 140 K, which was previously observed by 
/iSR (Ref. 1311 ) and neutron-diffraction^ 2 - measurements. 
In the doped NdFeAsO0.85F0.15 both the structural and 
AF transitions are absent at least down to the supercon- 
ducting transition temperature Tq = 43 K, below which 
the 75 As resonance rapidly disappears (Fig. [3]). Line- 
shape simulations also show that vq is almost temper- 
ature independent and therefore cannot account for the 
large temperature dependence of the observed 75 As NMR 
shift. 



IV. DISCUSSION 

In order to understand the observed 75 As resonance 
for the undoped NdFeAsO in the AF phase, we assume 
that this phase is characterized by the spin-density wave 
vector Q = (1,0,1) and that the ordered Fe magnetic 
moments are aligned along the a axis, i.e., similarly as in 
LaFeAsO 4 /iSR measurements suggested that the Fe or- 
dered moment in NdFeAsO is similar to that of LaFeAsO, 
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i.e., fiF e ~ 0.3 hb£^ The shift of the 75 As resonance in 
the AF phase is due to the internal field B mt aligned 
along the crystal c axis and is given by the anisotropic 
(off-diagonal) part of the hyperfine interaction with Fc 
moments, a xz , as B mt = Aa xz ^Fc>^ Its isotropic part in 
the ordered phase is filtered out from the shift due to the 
highly symmetrical position of 75 As with respect to Fc 
moments. The projection of this field adds to the external 
magnetic field, corrected for an isotropic hyperfine shift 
caused by Nd moments, B and gives rise to the effective 
magnetic field B e g = B Q + B lnt cost?. Here, $ denotes the 
angle between the external and internal magnetic fields. 
Taking into account also the first-, Ai/q\-d,(f>), and the 

second-order, Az/q <f>), quadrupole effects, the 75 As 
resonance frequency reads as 

v m ~m-i = 7A s 5cff + Ai/g> (0, 4) + Az/g } 0) . (2) 

Here, 7a s is the As gyromagnetic ratio and d and <f> 
define the orientation of B c g with respect to the prin- 
cipal axes of EFG tensor. Our simulation with Eq. 
(1) results in a good fit of the 80 K spectrum when 
o-xz = 10.5 kOe/fiB was used (inset to Fig. [5]). We stress 
at this point that the large anisotropic hyperfine fields 
cannot be explained by dipolar fields of the Fc moments 
and that they rather provide experimental evidence for 
the Fe 3d and As Ap hybridizations. 33 ' 34 
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FIG. 5: (Color online) Temperature dependence of the 75 As 
NMR shift measured for undoped NdFeAsO (open circles) 
and F-doped NdFeAsO0.85F0.15 powder samples (solid cir- 
cles). For comparison we show the corresponding shift mea- 
sured in LaFeAsOo.gFo.i (solid triangles) (Ref.ll7h. Inset: the 
75 As NMR line shape simulations [Eq. 1] for NdFeAsO sam- 
ple above (T = 180 K) and below (T = 80 K) the AF ordering 
transition Tat = 140 K. 

We now turn to the determination of the isotropic part 
of the 75 As hyperfine interaction measured by the Knight 
shift, K, in the paramagnetic phase (Fig. [5j . We stress 
that the measured shift is by an order of magnitude larger 
and has a completely different temperature dependence 



when compared to the one for LaFeAsOo.gFo.i M Tem- 
perature dependence of K in NdFeAsOi^Fj; can be well 
described by the Curie- Weiss-type behavior, K(T) = 
K + k/(T - 6). The value of the Curie- Weiss tem- 
perature 6 = -13(2) K [k = 0.99(2) K] is the same 
as determined from the magnetic susceptibility. We can 
therefore assume that the Curie- Weiss part of the shift 
comes from the Fermi contact hyperfine interaction be- 
tween the 75 As nuclei and Nd 4/ moments. The value of 
K = 0.096(15)% is close to the total NMR shift reported 
for LaFeAsOo.gFo.i^ We thus associate this part of the 
75 As NMR shift to the coupling of 75 As with the itiner- 
ant Fe 3d electrons in the FeAs layer. Interestingly, for 
undoped NdFeAsO we find that Kq is slightly reduced, 
i.e., K = 0.091(20)%. 

We proceed with the above analysis by recognizing that 
Kq is given by the hyperfine coupling to itinerant FeAs 
electrons and the orbital part, 

a iso 

A'o = -r^XFcAs + K olh . (3) 
Na^b 

Here Aorb = —0.075% is the orbital shift^i 7 - Na and /j,b 
are the Avogadro number and Bohr magneton, XFcAs is 
the molar susceptibility for the itinerant FeAs electrons, 
and finally a,p° As is the corresponding isotropic hyper- 
fine coupling constant. Since a precise dependence of 
XFeAs is in NdFeAsOi-^Fa; completely masked by the 
magnetism of the Nd 4/ moments (Fig. [5]) , we approx- 
imate it with the dependence for LaFeAsO taken from 
Ref. [I?]- Using the value of A' we derive ap° As = 
27.7 kOe/^ B and af° As = 29 kOe//i s for NdFeAsO and 
NdFcAsO0.85F0.15, respectively. Although these values 
have to be taken with certain caution due to approxima- 
tions used, we notice that they are about 10%-15% larger 
than in LaFeAsOo.gFo.i^ 

The contribution arising from the coupling to Nd 4/ 
moments, i^Nd, can be expressed as 

iso 

K Nd = K /(T-e) = -p^- X m, (4) 
NaUb 

where XNd = C/(T — 9) is the molar susceptibility of 
localized Nd 4/ moments and a^j is the corresponding 
isotropic hyperfine coupling constant. Using the value of 
k we calculate a}^ = 3.8 kOc//is. Although is by an 
order of magnitude smaller than ap° As , it is still surpris- 
ingly large as the REFcAsO system has so far been be- 
lieved to comprise well isolated Fe-As and REO layers^ 
The coexistence of itinerant Fe 3d and localized Nd 
4/ moments in NdFeAsOi- x F K is highly intriguing, and 
one wonders how the magnetic and/or superconducting 
properties are affected. Our experimental results clearly 
prove that the 75 As NMR probe interacts both with the 
Fe 3d conducting electrons in the FeAs layer as well as 
with the localized Nd 4/ moments in the NdO layer. The 
sizeable anisotropic hyperfine coupling constant between 
75 As and Fe 3d moments (a xz = 10.5 kOe/^s) in the AF 
phase has its origin in the hybridization between the Fe 
3d and the As Ap states. 
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FIG. 6: (Color online) Total density of states for LaFcAsO 
(shaded red) and NdFeAsO (full light blue line) using the ex- 
perimental room-temperature structural data. The density of 
states of the fictitious NdFeAsO) (labeled NdFeAsO La , using 
set 1 in Table. I) is shown by the dashed dark blue line. In- 
set: the density of states of the Fe 3d states for the same three 
structures. 



This result also nicely corroborates with the band- 
structure calculations and the density of states (DOS). 
In Fig. M the total DOS for LaFcAsO and NdFeAsO 
(using the room-temperature experimental structural 
parameters) and the fictitious Nd compound (labeled 
NdFcAsO La and calculated by using the experimental 
structure of the La compound at room temperature) is 
shown. The main features of the total DOS are very 
similar for the La and the fictitious Nd La compounds — 
except for the occupied Nd 4/ band (around —5 eV). 
Shorter lattice parameters in the Nd compound (go- 
ing from the dashed to the full line in Fig. [6|) have 
mainly an effect on the position of the 4/ and the Fe 
3d states (see below). Like in the La compound, we find 
also in the Nd compound a strong hybridization of Fe 
3d with As 4p states in the energy region between —4 
and —2 eV. The differences in these DOS can be seen 
more clearly on a smaller scale. The Fe 3d states are 
slightly shifted in direction of lower energy from the La 
via the fictitious Nd La to the Nd compound; see inset 
of Fig. [6j These subtle electronic structure changes are 
also observed in the 75 As NMR experiment. We no- 
tice that the As-Fe isotropic hypcrfinc coupling constant 
af? c ° As = 27.7 to 29 kOc/^ B is about 10% larger than that 
reported for the La-family. This constant depends on 
the degree of polarization of the As 4s electrons caused 
by the As 4p states via exchange polarization effect ^ 
Therefore it may be — similar to the DOS — sensitive to 
the small crystalline structure variations such as lattice 
contraction or the shift of the As along the z axis found 
for NdFeAsOi-^Fa; in our x-ray diffraction (XRD) ex- 
periments. We conclude that all 75 As NMR parameters, 



a IZ , a,p° As and vq arc very sensitive to the details of As 
4p electron distribution. Related to that, we refer here 
to the work of Mukuda et ali^i where the monotonic in- 
crease in the superconducting transition temperature Tc 
with 75 As quadrupole frequency vq has been reported. 
Based on these observations, our calculations (Table I) 
may suggest that it is the lattice contraction and not 
the presence of additional RE 4/ states that shifts Tc to 
higher values in NdFcAsOi_ x F x compounds. This is also 
in line with pressure experiments on LaFeAsO, where T c 
was increased up to 43 

Finally we comment on the non-negligible coupling be- 
tween 75 As and Nd 4/ moments, a^° A = 3.8kOe//j,s, 
which directly evidence an interaction between the REO 
and FeAs layers. Such interlayer coupling could open the 
possibility for an indirect exchange interaction between 
the localized Nd 4f moments mediated by the Fe 3d 
and/or Nd 5d (partially occupied in the band-structure 
calculations) conducting electrons, i.e., a Ruderman- 
Kittel-Kasuya-Yosida (RKKY)-type interaction. The 
small Curie- Weiss temperature (8 = —14 K) as well as 
the low ordering temperature for the Nd moments may 
be due to the long-range nature of the RKKY-type in- 
teraction and its oscillatory dependence on the distance 
between the Nd moments. Additional experiments are 
needed to shed more light on the size of Nd ordered mo- 
ment at low temperatures. 



V. SUMMARY 

In summary, we have performed a study of the struc- 
tural, electronic and magnetic properties of the super- 
conducting NdFeAsOi-^F* (T c = 43 K for x = 0.15) 
by high-resolution synchrotron x-ray powder diffraction, 
magnetization, and 75 As NMR measurements. Struc- 
tural investigations showed a contraction of the lattice 
structure with shorter RE-As bond distances, an in- 
creased thickness of the As-Fe2-As layer and less dis- 
torted Fc-As4 tetrahedra for the NdFeAsO family when 
compared to the LaFeAsO family. A typical 75 As 
quadrupole powder line shape measured for NdFeAsO 
yielded an enhanced quadrupole frequency with respect 
to the La analogs. Furthermore, we calculated the elec- 
tric field gradient using the band-structure code FPLO 
and obtained good agreement for the undoped com- 
pound. The enhanced EFG is due to a more prolate 
As 4p electron distribution, which is mainly caused by 
reduced lattice parameters and not the presence of Nd 
4/ electrons. We also found a scaling of the 75 As NMR 
shift with the magnetic susceptibility of localized Nd mo- 
ments suggesting a weak coupling between the NdO and 
FcAs layer. The RKKY- interaction between localized 
Nd 4/ moments mediated by itinerant Fe 3d and/or Nd 
5d states may originate from such a coupling. 
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